Thermodynamic properties of Pb determined from pressure-dependent critical-field 

measurements 
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We have carried out extensive low-temperature (1.5 to 10 K) measurements of the critical field, H^, 
for the element Pb up to a pressure of P = 1.2 GPa. From this data the electronic entropy, specific 
heat, thermal expansion coefBcient and compressibility is calculated as a function of temperature, 
pressure and magnetic field. The zero-field data is consistent with direct thermodynamic measure- 
ments and the P-dependence of Tc and specific heat coefficient, 'y{T, P) allows the determination of 
the P-dependence of the pairing interaction. 

PACS numbers: 74.25.Bt, 74.62. Fj, 74.25.Ha, 74.25.Ld 



The thermodynamic and superconducting properties 
of the metal Pb have been extensively studied, in- 
cluding specific heatii^, isotope effect^i^, critical field^ 
and thermal expansion®. In addition the effects of 
pressure on its superconducting properties have also 
been well studied^i^ii^. Subsequently, the use of the 
pressure-dependent superconducting transition temper- 
ature Tc{P) for Pb as a low-temperature manometer 
has also been propose d^°i^^'^^ . The above quoted ref- 
erences are only the early ones and one would assume 
that the thermodynamic properties of this material have 
been comprehensively investigated. 

In the course of investigations into the pressure depen- 
dence of the oxygen isotope effect on Tc for the high-Tc 
superconductor YBa2Cu408 we used the metal Pb as 
an internal pressure calibrant. Using a clamp cell in a 
SQUID magnetometer we tracked the pressure using the 
reported pressure dependence of Tc for Pb--. Despite the 
many above-noted early studies on Pb we discovered that 
the pressure, P, and temperature dependence of the crit- 
ical field. He, for Pb has not been reported except at very 
low pressures (< 0.03 GPa)^. We therefore carried out a 
study of Hc{T,P) from 1.5 to 10 K and up to 1.2 GPa 
from which we have derived a full set of thermodynamic 
parameters. This paper presents a summary of these 
measurements and the deduced values of the T- and P- 
dependence of the electronic entropy, S{T,P), electronic 
specific heat coefficient, ^{T,P), thermal expansion co- 
efficient, a{T,P) and electronic compressibility, k{T,P). 

The sample used was a bar of dimensions 1.1 x 
0.93mm cross section and 5.58mm long, shaped from Pu- 
ratronic 99.9985%, 2 mm diameter Pb wire from Johnson 
Matthey. Using the formulae of Osborn the demagnetiza- 
tion factor was calculated as 0.97 and this value was used 
throughout. The Pb sample was annealed in vacuum at 
300 C for 8 hours to reduce magnetic hysteresis^. This 
was loaded in a 2.67 mm diameter 9 mm long Tefion cap- 
sule along with Fluorinert FC70 and FC77 mixed in 1:1 
ratio as a cryogenic hydrostatic pressure medium. The 
sample capsule was then placed in a miniature home- 
built non-magnetic Be-Cu(Mico Metal 97.75% Cu, 2% 
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FIG. 1: (Color) Magnetization curves for the Pb sample plot- 
ted for various applied fields (annotated in Gauss as 10, 100, 
200 etc) and four different pressures (0, 0.5, 0.9 and 1.15 GPa). 



Be) piston clamp cell (8.8mm diameter, 65mm length, 
with cobalt-free tungsten-carbide pistons (Boride). The 
pistons of this cell are lightly tapered using electric- 
discharge machiningi^. To apply pressure the cell was 
preloaded before clamping at room temperature using a 
laboratory press with calibrated digital pressure gauge 
(Ashcroft Model 2089, 0.05% accuracy). The magnetiza- 
tion measurements were carried out in a Quantum Design 
MPMS SQUID magnetometer. The pressure in the sam- 
ple was measured from the reported shift in Tc of Pb 
at zero fieldi^. For in-field measurements the cell was 
always zero-field cooled to avoid hysteresis error. 

Fig. 1 shows the measured T-dependence of the mag- 
netic moment, M, up to a field of 700 G and at pressures 
of GPa (1 bar), 0.5, 0.9 and 1.15 GPa. Tc was deter- 
mined by the extrapolation to zero of the steepest slope 
of M{T). This gives Tc{H) which is rc-plottcd in Fig. 
2 as Hc{T,P) versus T. The data shows a progressive 
decline in both Tc and Hc{T) with increasing pressure. 
We will see that the combination of these two pressure 
effects allows access to the pressure dependence of the 
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FIG. 2: (Color) The T-dependence of the critical field, He 
for pressures of 0, 0.5, 0.9 and 1.15 GPa. Data points are 
measured data from Fig. 1, curves are fits using eq. 2. 

density of states, N{Ep), at the Fermi level, Ep and of 
the pairing interaction. 

To fit this data we expand Hc{T) as a power series in 
T . The first-order term must be zero otherwise 5^ 7^ 
at T = and the third-order term must also be zero to 
avoid a negative 7 at low T. It is usual, therefore, to 
adopt an even polynomial of the forni^ 

H,{T, P) - H,o{P) + H,i{P) xT^ + H,2{P) X (1) 

where, as indicated, the coefficients are each an indepen- 
dent function of pressure. As we will see. He cannot be 
linear in T because of the third law requirement that the 
entropy, S* ^ as T — > 0. It turns out that the coef- 
ficients are not independent and it is common to adopt 
the "similitude principle" which separates the P- and T- 
dependent terms as follows^ 

i?e(T,P) = H,oiP) X fit) 

= Hcq{P) X {l-axf + {l-a) X t^}{2) 

where t — T/Tc- A free fit of eq. (1) to the data consis- 
tently yielded Hc2/Hco very close in value to l—Hd/Hco, 
consistent with eq. (2). We therefore proceeded to fit the 
data with eq. (2) for all subsequent analysis. The fits are 
shown by the solid curves through the data shown in Fig. 
2. We obtain Hco values of 804.08, 770.88, 739.99 and 
729.01 Gauss for P = 0, 0.5, 0.9 and 1.15 GPa, respec- 
tively. A linear fit yields 

HcoiP) = (803.72 - 67.028P) Gauss (3) 

This now allows the scaling of all the data by plotting 
Hc{T,P)/Hco{P) versus as shown in Fig. 3. The 
small departure from linearity in Fig. 3 reflects the small 
quartic term (a = 0.95475; 1 - a = 0.04525) and the 
scaled data in Fig. 3 allows a global fit to determine 
a that provides an overall consistency in the data fits 
that is necessary when second derivatives of Hc{T, P) are 
used to determine the electronic specific heat coefficient, 
thermal expansion coefficient or compressibility. 
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FIG. 3: (Color) Scaling of the T-dependent critical field, 
showing Hc{T, P) /Hc{T = 0, P) versus t'^ = (T/Tcf for pres- 
sures of 0, 0.5, 0.9 and 1.15 GPa. The solid curve is the func- 
tion /(T) = 1 - 0.95475t^ - 0.04525**. Inset: 7„ calculated 
from eqn. (6). 

These parameters may be determined as follows. We 
restrict our thermodynamic parameters to the electronic 
contribution, thus ignoring the lattice contribution to the 
specific heat or thermal expansion. The difference in free 
energy between the normal and superconducting states 
at the same T and P is given by 

G,(T, P) = G„(T, P) - ^^oH,{T, Pf X Vm (4) 

where /io is the permeability of free space and Vm is the 
molar volume. As shown by Shoenbergi^ the thermo- 
dynamic parameters are obtained by differentiation with 
respect to T and P, giving 

Ss[T,P) = S,,iT,P) + fioHcVM{dHjdT)p 

Is = Jn+t^oVMidH,/dT)l 

- ^ioHcVM{^^HJ^T^)p 

V, = Vn-floHcVM{dHjdP)T 

as = a„-fioidHjdT)px{dHjdP)T 

- ^iM^^'HJ^T^P) 

Ks = Kn+ fJ-o{dHc/dP)p 

+ ^ioH,id^H,/dP^)T (5) 

Here terms such as {l/2)^oH^VMCttot in the entropy or 
HoHtotHc{dHc/dP)T in the compressibihty obtained by 
differentiation of the molar volume are ignored as negli- 
gible, where atot is the total volume coefficient of ther- 
mal expansion including lattice and electronic terms, and 
similarly for Ktot- 

Taking the first of these expressions, that for the en- 
tropy Ss(T,P), and dividing through by T we note the 
requirement that Ss/T — > as T — > due to the opening 
of a full superconducting gap in the DOS. At the same 
time if 7„ is constant then Sn/T is just 7„. This then im- 
poses a relationship between 7„ and the quadratic term 
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FIG. 4: (Color) The T-dependence of (a) the electronic en- 
tropy in the normal state and superconducting state for pres- 
sures of 0, 0.5, 0.9 and 1.15 GPa; and (b) of S/T. 



in Hc{T). In particular, 



7„ = lim -fioT~^H,VMidH,/dT)p 



(6) 



In this way 7^ is calculated from the fits in Fig. 3 
and its value is plotted in the inset to Fig. 3. The value 
at ambient pressure, 7„ = 3.47 mJ/(mol K^) compares 
favorably with the value 7„ — 3.13 mJ/(mol K^) de- 
termined from direct specific heat measurementsi. The 
pressure dependence of 7„ shown in the inset reveals a di- 
mensionless volume dependence given by d\njn/d\iiV = 
3.26 ± 0.27 (where we assume a total compressibility of^^ 
kq = 0.0205 GPa^^). There have been various estimates 
of this dimensionless parameter ranging from 1.7 ± 0.5 
from thermal expansion measurements^^, 1.8 ± 0.5 from 
volume expansion at T^ii, 6.0 based on pressure depen- 
dence of and 3.1 ± 0.8 in good agreement with our 
present value based on more recent measurements of vol- 
ume change at Tc by Ott^^. Note that this is consid- 
erably stronger than the value for the free electron gas 
d\njn/dlnV = 2/3. 

With these relations in place we proceed to calculate 
the entropy as shown in Fig. 4(a). The linear slope in the 
normal state is just the above-determined 7„ (P) . The en- 
tropy in both states must vanish at T = 0. A more rigor- 
ous test of the polynomial fits to the data is whether Sg/T 
also vanishes as T ^ and whether the data preserves 
monotonic systematics in this region. The T-dependence 
of S/T is plotted in Fig. 4(b) and this indeed shows a 
quadratic behavior at low T which is perfectly system- 
atic with increasing pressure and extrapolates to zero as 
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FIG. 5: (Color) The T-dependence of the specific heat coef- 
ficient 7(r, P) in the normal and superconducting states for 
pressures of 0, 0.5, 0.9 and 1.15 GPa. The small departure 
from monotonic systematics in the specific heat jump for 0.9 
GPa reflects the accumulation of errors in the 2"'' derivatives. 



T — * 0. In a finite external field the same entropy curves 
are retraced up to the reduced Tc{H) value and then the 
entropy jumps discontinuously to the normal-state value, 
consistent with a first-order phase transition in magnetic 
field and second-order when H — 0. 

The condensation energy, Uo is determined from 



[Sn{T,P)~Ss{T,P)]dT 



(7) 



from which we obtain the following values: Uo — 47, 
43.2, 39.8 nd 38.6 mJ/mol for P = 0, 0.5, 0.9 and 1.15 
GPa, respectively. The condensation energy reduces with 
pressure due to the twin effects of reduced Tc and reduced 

In- 

Turning to the specific heat coefficient, this is calcu- 
lated using the second equation in eqs. (4) and is shown 
in Fig. 5 for both the normal and superconducting states. 
The entropy balance is evident from the equal areas 
above and below the temperature where 7^ = 7n, and is 
confirmed by the fact shown in Fig. 4(b) that SsjT — 7„ 
at Tc. We find that the jump in 7(T, P) at Tc is A7C = 
7.85, 7.51, 7.05 and 7.17 mJ/(mol K^) for P = 0, 0.5, 0.9 
and 1.15 GPa, respectively. Values of A7c/7n are respec- 
tively 2.33, 2.28, 2.18 and 2.24. There is a weak pressure 
dependence here in this ratio that is not expected with a 
simple BCS picture where A7c/7„ = 1.55. Direct specific 
heat measurements at ambient pressure have variously 
obtained A7C = 7.33 ±0.3 mJ/(mol K^) (Clement and 
QuinneliS), 7.96 ±0.08 mJ/(mol K^) (Shiffman et al.^") 
and 8.13 ±0.13 mJ/(mol K^) (Neighbor et al.^^) in good 
agreement with our ambient result. Combining with the 
measurement by Horowitz et ali- of 7^ = 3.13 mJ/(mol 
K^) gives A7c/7„ — 2.34, 2.54 and 2.60, again in good 
agreement with our present ambient-pressure result. We 
note that for the purely parabolic form He — H°{1 — t^} 
then A7c/7n = 2 exactly, and so the experimentally ob- 
served excess of A^/c/jn is a direct indication of the pres- 



4 



CD 

o 

I 



o 
E 

E 
o 



0.2 



0.1 - 



0.0 



-10 



-0 GPa 
-0.5 GPa 
0.9 GPa 
-1.15 GPa 




012345678 
T [K] 

FIG. 6: (Color) (a) The T-dependence of the difference in 
volume thermal expansion coefficient between the supercon- 
ducting and normal states for pressures of 0, 0.5, 0.9 and 1.15 
GPa. (b) the difference in molar volume for the same pres- 
sures. 



ence of a higher order (quartic) term in the T-dependent 
critical field. 

Finally, Fig. 6(a) shows the calculated difference in 
volume thermal expansion coefficient between the super- 
conducting and normal states at each pressure and Fig. 
6(b) shows the difference in molar volume. Given that 
the absolute molar volume at T = K is 18.26 cm'^/mol 
these changes in volume seem extremely small but they 
are readily measureable. The jumps in thermal expansion 
coefficient at Tc in zero field are Aa = 1.91, 1.63, 1.42 
and 1.29 x 10"^ for P = 0, 0.5, 0.9 and 1.15 GPa, 
respectively. We omit the calculations of the electronic 
compressibility because these involve double derivatives 
with respect to pressure and here the errors begin to ac- 
cumulate markedly. Nonetheless the discontinuous jumps 
in isothermal compressibility may be calculated from the 
Ehrenfest equation for a second order transition: 



dTc/dP = Ak/Aq; 



(8) 



Empirically we find that 

TciP) = 7.1863 - 0.3847P + 0.01769p2 (9) 

so that the absolute jumps in isothermal compressibil- 
ity at Tc in zero field are Ak — -7.3, -6.0, -5.0 and -4.4 
X 10"^ GPa-i for P = 0, 0.5, 0.9 and 1.15 GPa, re- 
spectively, while the relative jumps are Ak/kq — -3.56, 
-2.93, -2.44 and -2.15 ppm. Direct measurements^^ of the 
ambient-pressure elastic moduli for Pb in the supercon- 



ducting state yield Ak/kq = -4.0 ppm, again in excellent 
agreement with our analysis for P = GPa. 

We conclude by using the above results to estimate 
the pressure-dependence of key parameters in the BCS 
model^^. The condensation energy is 



-NiO)NAAl 

(47.11 - 7.75P) mJ/mol (10) 



where P is in GPa. The jump in specific heat is given by 

Is - In 



10.2fc|7V(0)iVA 
(7.85 - 0.68P) mJ/(mol K^). (11) 



Therefore, we obtain 

N{0) = (1.07- 0.093P) states/atom/eV(12) 
Ao = (11.1 -0.432P) K (13) 
2Aa/kBTc ^ 3.09 + 0.037P. (14) 

This should be compared with the BCS result 
2Ao/kBTc = 3.53 and the value of Aq = 15.76 ± 0.05 
K observed from tunnelingSi giving 2Ao/fcB Tc = 4.37. 

If, on the other hand, we determine the DOS from the 
nearly-free electron expressioiii^^ 

7„ = i2/3)7:^klNiO)NA 

= (3.47 - 0.232P) mJ/(mol/K2). (15) 

Then we obtain 

iV(0) = (0.736 - 0.0492P) states/atom/eX16) 
Ao = (13.7 - 0.653P) K (17) 
2Ao/kBTc = 3.72 + 0.007P. (18) 

and the ratio 2Ao/fcBTc is more or less pressure indepen- 
dent. Finally, using the BCS expression^ for Tc 



Tc = 0.85613 cxp(- 



1 



N{0)W 



(19) 



where Qd is the Debye temperature (= 96 K for Pb2^) 
and W is the pairing interaction, then, using eq. (12) 



W = (0.343 + 0.0167P) eV 
or using eq. (16) 

W = (0.499 + 0.0143P) eV. 



(20) 

(21) 
-2.6 



Here we have used the Gruneisen coefficient 
for Pb25.. 

In summary, we have measured the pressure depen- 
dence of the superconducting critical field in Pb to 1.15 
GPa using a clamp cell in a SQUID magnetometer. By 
using thermodynamic identities we have determined the 
electronic free energy, entropy, specific heat coefficient, 
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thermal expansion coefficient and compressibiHty, includ- 
ing the jumps in these properties at T^.. The calcu- 
lated results match rather well the experimentally ob- 
served ambient pressure data, where it is available. The 
present results allow calculation of the these parameters 
as pressure-dependent quantities and an estimation of the 



pressure dependence of the density of states, supercon- 
ducting energy gap and the pairing interaction. 
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